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Instanton Distribution in Quenched and Full QCD 
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In order to optimize cooling as a technique to study the instanton content of the QCD vacuum, we have 
studied the effects of alternative algorithms, improved actions and boundary conditions on the evolution of single 
instantons and instanton anti-instanton pairs. Using these resuhs, we have extracted and compared the instanton 
content of quenched and full QCD. 
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1. Introduction 

Because of the significant role that instantons 
play in light hadron structure |^J^ and the in- 
trinsic importance of understanding topological 
excitations in QCD, it is of fundamental interest 
to understand the instanton content of the QCD 
vacuum and, in particular the effect of dynamical 
quarks on it. Although cooling is a powerful tech- 
nique to preferentially filter out non-topological 
excitations relative to instantons, it is limited by 
the removal of instantons due to lattice artifacts 
and instanton anti-instanton annihilation. We 
address this problem and compare quenched and 
unquenched results. 

2. Relaxation Algorithm 

In order to extract the instanton content of the 
lattice configurations efficiently on a parallel com- 
puter, we used a variant of the cooling method by 
discretizing the relaxation equation. 
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introducing a small step size parameter e in the 
"relaxation time" r and updating simultaneously 
all links on the lattice. For large values of e this 
algorithm is unstable but in the limit e — > it 
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converges to the solution of the relaxation equa- 
tion(0). We have found that for SU(3) gauge field 
the value e = 0.025 gives fast and stable relax- 
ation and we have used this value in our measure- 
ments. Comparing the evolution of the action for 
our relaxation method with e = 0.025 and the 
"standard" Cabbibo-Marinari cooling, we found 
that one cooling step is approximately equivalent 
to 4 relaxation steps and that the cooling histories 
are very similar. 

3. Identifying Instantons 

We used the following procedure to identify in- 
stantons in our configurations. We consider all 
peaks in the action and topological charge density 
as candidates for the instantons. With an initial 
value of the instanton position xq (peak position) 
and instanton size po — (48/5*0)^/^ (with the nor- 
malization S{x) — s- a^FF in the continuum), we 
perform a least squares fit to the action and topo- 
logical charge densities of a classical continuum 
instanton. 
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as a function of xo,p. If the fitting process con- 
verges and the final values for the instanton posi- 
tion and size p are within a pre-defined range 
close to the initial values, we record the identi- 
fied instanton. If one of the conditions above is 
not satisfied, the candidate peak in the action or 
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topological charge is discarded. 

Full details of the algorithm will be published 
elsewhere Although there is some degree of 
arbitrariness in selecting the parameters in the 
algorithm, we were able to identify around 50% of 
the peaks in the action distribution as instantons 
with the remaining peaks corresponding to small 
instantons or other excitations. 

4. Isolated Instantons 

Several effects can distort the instanton dis- 
tribution during the relaxation process. The 
smallest instantons may disappear by falling 
through the mesh. The larger instantons may 
shrink because of interaction with periodic images 
or may disappear by instanton anti-instanton 
pairs annihilate. We separated these effects 
by studying the evolution of discretized semi- 
classical instanton configurations. 

First, we have investigated the shrinking of a 
single isolated instanton for a variety of instan- 
ton sizes and two lattice sizes, 16^ and 24"*. We 
studied two actions, the Wilson action and a first 
order improved action [Q, 

^imp = jTr(l - Wixi) - ^Tr(l - ^2x2). (3) 

The first order term in the Wilson action expan- 
sion in I is negative and the relaxation pro- 
cess forces instantons to shrink. The coefficients 
in the improved action were chosen so that the 
first order term is zero and the second order term 
is positive. In Fig. |l|, we show the evolution of a 
single isolated instantons with sizes p = 3.0 and 
p = 7.0. Notice that during relaxation the trajec- 
tories terminate for the smaller instantons as they 
"fall through the lattice" and that the small in- 
stantons are much more stable with the improved 
action jj] . Adding more terms to the action could 
further stabilize instantons ||] but the first-order 
improved action was simple and already sufficient 
for our studies. For large instantons the evolution 
does not depend appreciably on the choice of the 
action but the boundary effects are substantial. 
Since we are ultimately interested on the order 
of 100 relaxation steps at most, we conclude that 
with the improved action and significantly larger 



lattices, the artificial loss of single instantons in 
cooling is negligible. 
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Figure 1. Evolution of isolated instantons size 
p = 3.0 and p = 7.0 under relaxation for Wilson 
and improved actions on lattice sizes 16* and 24*. 

Next, we studied the evolution of an instan- 
ton anti-instanton (I-A) pair. Figure || shows 
the history in a typical case: instantons of ini- 
tial sizes Pi = PA = 6.0a with the separation 
between centers s = 12.0a. On the small lat- 
tice 16*, the boundary effects are much stronger 
than the interaction between objects and they 
shrink in place before they have a chance to inter- 
act. On the larger lattices the instanton and anti- 
instanton attract each other and annihilate. The 
finite volume effects have stabilized on the largest 
lattices. We have also checked that the discretiza- 
tion effects are negligible by comparing the evolu- 
tion of the two appropriately rescaled configura- 
tions. (In Fig ^ the "data" for the 36* lattice are 
actually obtained by a re-scaling from a 24* lat- 
tice.) 

The "interaction" time during which the in- 
stanton move toward each other is much less then 
the "shrinking" time and the individual instan- 
tons in a pair do not shrink significantly. In the 
case of a real dynamical configurations, we ex- 
pect the instantons to interact more strongly and 
hence we believe the conventional use of periodic 
boundary conditions will not introduce serious fi- 
nite volume distortions. 
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Figure 2. Evolution of I- A pair under relax- 
ation on the lattices of different size. Upper 
curves gives the I-A separation s and the lower 
the sizes p. 



5. QCD Results 

We have studied 21 full QCD configurations 
for (3 = 5.5, 2 flavors of Wilson fermions and 
K = 0.160. As measured by LANL group dj, 
this choice of parameters corresponds to the lat- 
tice spacing a(/^) = O.llfm and pion mass — 
360MeV. Configurations were separated by 50 
HMC trajectories with average length 50 steps 
of size e — 0.01. Those configurations have the 
same value of Kc = 0.16145 as quenched QCD for 
13 = 5.858 (See Ref. @). To compare quenched 
and full QCD in roughly the same physical region, 
we have generated 23 quenched configurations for 
(3 = 5.85 separated by 500 heat-bath iterations. 
We checked also that the Creutz ratios in both 
cases are comparable in the range of 3-5 lattice 
units. 

Figure |^ shows the distribution of instantons in 
sizes after 20, 30 and 50 relaxation steps. Since 
our method for identifying instantons is not re- 
liable before 20 steps, we can not determine the 
early evolutions of the distribution, but after 20 
the distribution certainly changes dramatically 
with relaxation. Given that the number of instan- 
tons are 0(200) or larger in a 16* box, the aver- 
age distance between instantons is about 4 lattice 
units, comparable to the average instanton size. 
Under these conditions the interaction between 



instantons is very strong and consequently the 
effects of TA pair annihilation are much stronger 
than boundary effects. 

Since the erosion of the instantons distribution 
with cooling should be the same for the quenched 
and unquenched samples, we believe the fact that 
the distributions shown in Fig.^ are essentially 
identical within errors provides strong evidence 
that the physical instanton distributions are very 
similar in quenched and full QCD, at least at this 
sea quark mass, which is of the order of m^. 
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Figure 3. The instanton distribution in quenched 
and full QCD. The smooth curve is the measured 
distribution convoluted with a Gaussian curve of 
width O.Olfm. Error bars are estimated by the 
jackknife method. 
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